Optically induced orientational transitions in nematics with planar alignment 
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Theoretical study of dynamical phenomena induced by a linearly polarized plane wave incident 
perpendicularly on a planar aligned nematic layer with the light intensity as the control parameter 
is reported. We find the threshold of the Optically Induced Freedericksz Transition for the planar 
state as a function of the problem parameters. It occurs that the threshold is substantially lower 
than that expected before. Above the primary instability the director settles either to a stationary 
or to an oscillatory states depending on a thickness of the layer. These states become unstable at a 
secondary threshold through a heteroclinic bifurcation and the director settles to a new stationary 
distorted state. 

PACS numbers: 42.70.Df, 05.45.-a, 42.65.Sf 



Optical phenomena exhibited by nematic liquid crys- 
tals has been a subject of intensive study during the last 
few decades. A nematic behaves optically as a uniaxial 
anisotropic medium with the optical axis along the local 
molecular orientation described by the director n{r,t). 
The dielectric tensor governing the propagation of light 
is anisotropic and depends on n. There are two compet- 
ing mechanisms which determine the nematic alignment. 
On one hand, it is enforced by the anchoring forces at 
the confining surfaces. On the other hand the electric 
field of the incident light exerts a torque which may con- 
flict with the boundary-imposed alignment. The torque 
increases with an increase in the light intensity. Even- 
tually, at a certain light intensity the boundary-imposed 
alignment becomes unstable being replaced by the light- 
imposed one (Optically Induced Freedericksz Transition) 
which in turn affects the light propagation A large 
variety of nonlinear phenomena occurs as a result this 
feedback 0. The possibility of dynamical changes of 
the nematic refractive index employing the orientational 
phenomena, has attracted recently much attention in the 
context of construction of all-optical devices based on 
photonic structures infiltrated with liquid crystals, see 
e.g. and references therein. The goal of this paper is 
elaborated theoretical study of light induced phenomena 
in a nematic cell with planar alignment (without pho- 
tonic structure) which is much less explored than that 
with the homeotropic one. This problem was considered 
before in [1] where a linear analysis of the planar state 
with respect to a restricted class of perturbations was 
done. In the present paper we perform the general linear 
analysis with respect to arbitrary (small) perturbations. 
We show that the planar state might lose stability either 
through a stationary bifurcation or through a Hopf bifur- 
cation when the thickness of the layer is above a certain 
critical value. The thresholds for the primary instabil- 
ity are substantially lower than those predicted before. 
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FIG. 1: Geometry of the setup: linearly polarized light along 
the y direction incident perpendicularly on a nematic layer 
with the director no || x (planar state). The components of 
the director n are described in terms of the twist angle <1> 
(<I> = in the planar state). 

Above the Hopf bifurcation the system settles to an os- 
cillatory state. This state, in turn, becomes unstable at 
some higher value of intensity where the transition to a 
new stationary distorted state via a heteroclinic bifurca- 
tion occurs. 

We consider a linearly polarized along the y direction 
plane wave, propagating along the z axis. The wave in- 
cidents perpendicularly to a nematic layer, see Fig. [T] 
Initially the director is orientated along the x axis (pla- 
nar alignment) and keeps its orientation at the bound- 
aries (rigid boundary conditions). We restrict the con- 
sideration analyzing only spatially-homogeneous in the 
layer plane solutions which depend on the z coordinate 
solely. To describe the director orientation the twist an- 
gle $(z,t) is introduced, so that n = (cos $, sin<I>, 0). 
The starting point is the set of nematodynamic equa- 
tions coupled with Maxwell's equations for the propaga- 
tion of light [H]. The dynamical equation of motion for 
^{z,t) is obtained from the balance of torques (elastic, 
electromagnetic and viscous) acting on the nematic. The 
components of the electric field might be represented in 
terms of the amplitudes of the ordinary and extraordi- 
nary waves Ae , that vary slowly with z on the scale 
{kt)L)~^ (see e.g. [Qj), where /cq — i^/c is the wavenumber 
of the incident wave and L stands for the layer thickness. 
Finally, in the absence of a velocity field the coupled PDE 
for (f> and ordinary differential equations (ODEs) for Ag 
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and Ao can be written as follows 



(1) 



= -(a,$)e'^-«Me, a.Ae = (a,$)e-*^-»M<„ (2) 

where fco is dimensionless thickness fco = 2L5n/X with 
(5ri = rie — rio ('T-e and n,, are, the refractive indices of 
the ordinary and extraordinary light, respectively) and A 
is the incident light wavelength. To make the variables 
in Eqs. H]), ^ dimensionless they have undergone the 
following scale transformation: t t/r; z ttz/L; am- 
plitudes Ae^o have been normalized over the electric field 
amplitude of the incident light; p — I / Ic, where / stands 
for the incident light intensity. Here r = ^iLp' / tt"^ K2 is 
characteristic relaxation time of the director motion and 
Ic = 87^^ci^'2^^e'5r^/A^(r^e + Uo), K2 is the twist elastic 
constant and 71 = 03 — 02 is the rotational viscosity. 

The boundary condition for $ and initial conditions 
for Ao, Ae at z = read as follows: 



*.=o,.W=0, |A„o|' = l, |AeoP = 0, A 



eO 



A* 



0. (3) 



Note that owing to the reflection symmetry with respect 
to the y direction Eqs. ([TJlS]) are invariant under the 
transformation S:{^, A^, Ao} {— $, =F^e, ±^0} • It is 
convenient to represent the solution of Eqs. ((I])-® as 
a series ^{z,t) ~ ^n{t) sin(nz), where each term of 
the sum satisfies the boundary conditions Eq. ^ identi- 
cally. Then, Eqs. (H])-© are transformed into an infinite 
set of coupled equations for To make the problem 
tractable the set is truncated at a certain large enough 
number of equations N which are solved numerically by 
the standard Runge-Kutta method. The error caused by 
the truncation is controlled by test runs with double and 
triple number of the modes. For every set of parameters 
N is selected so that the difference between the routine 
and test runs is better than 1%. Regarding equations for 
Ao,Ae, they have to be solved dynamically at each step of 
numerical integration for time t. In addition, we perform 
a linear analysis (numerically) of the stationary distorted 
states {dt^n = 0) by calculating the corresponding Jaco- 
bians. 

The starting point of the study is the linear stability 
analysis of the planar state ($(z,i) = 0). The linearized 
integro-differential equation for ^{z,t) — $(z)exp((Tt) is 
as follows: 

dl<^ + 2pkl (^* + ^o^ $(z')sin[fco(2' - z)]dz'j =o-$.(4) 

It results in the following set of equations for $„ (eigen- 
value problem): 



(5) 
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1 n^'5mn + 



pk^mn sin(7r/co) 
7r(m2 — fcQ)(n2 — fcg) 
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FIG. 2: Re(cr) and Im{a) vs. p for ko = 1.6. Filled cir- 
cle: primary threshold {pi = 0.95, Hopf bifurcation). Shaded 
tongues cross abscissa at, respectively, pi = 16.6, P2 = 17.8, 
pi — 24.3 and pi = 26.5 depicted by empty circles. 

The solvability condition for Eq. ([5]) requires det(A — a I) 
— 0. It brings about infinite number of eigenvalues (T„. 
However, at n » fcg or (and) ItQ the off diagonal el- 
ements of matrix Amn decay as 1/n or (and) l/m, while 
the diagonal elements equal — n^. Then, in the leading 
approximation det(A — a I) factorizes for a product of in- 
finite number of the diagonal elements at n > (which 
give stable real eigenvalues (T„ w — related to elastic 
relaxation) and a minor N x N, where N ^ k^. Eigen- 
values of the minor are related to light-induced pertur- 
bations and should be inspected more carefully. These 
arguments provides us with the natural truncation scale 
~ fcp. Note that the minimal number of modes for ad- 
equate description of the phenomenon increases as L^. 
The planar state looses stability when the real part of (at 
least) one of the eigenvalues i?e((T„) becomes positive. It 
should be stressed that matrix Amn is not Hermitian and 
may have complex eigenvalues, in other words unstable 
modes may be oscillatory. 

Typical results of the discussed stability analysis are 
presented on Fig. [21 where i?e(cr„) and Im{an) versus 
p are shown for kg ~ 1.6. (The calculation are made 
for N = 20.) It is seen that i?e(cr„) forms a tongued 
structure with a family of branches which, at first, appear 
as pairs of purely real branches that go almost parallel 
to each other as p increases, but then the pair merges at 

(n) 

a certain p>p\op- Next, the branches continue as a pair 
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of complex conjugate modes. Some of the purely real 
branches cross the zero line at intensities p given by the 
formula 

2p- sin[7rfcoVl + 2p\ ^ 
7rfco\/l + 2/0 

which was derived in Q . The fist two tongues depicted by 
shaded areas crosses the abscissa at, p\ = 16.6, P2 — 17.8, 
p\ — 24.3 and p\ = 26.5, respectively, and are the same 
as ones shown in Fig. 1 of (The dimensionless X = 5 
introduced there corresponds to ko — L/i: ^ l.Q used in 
our calculations.) In Q such a structure was interpreted 
as a series of alternate stable and unstable intervals for 
the planar state as the light intensity increases with the 
lowest threshold for instability given by p\ = 16.6. How- 
ever, in reality, the planar state looses stability via a Hopf 
bifurcation at much lower value pi — 0.95 (filled circle on 
Fig. [2|) and never restores stability at higher intensities 
again. In fact, the reason for such a difference with the 
results of [4?| is that for solving the eigenvalue problem 
(5) therein a particular (not a general) form of a trial 
solution was used [see [3, Eqs. (6,10)]. As a result, the 
growth rates turned out to be a priory real instead of 
being complex and, thus, the ones with Im{a) ^ pre- 
dicted here have been not caught. The stability diagram 
in the (fco , p) plane is shown in Fig. [3] by solid lines and 
the thresholds calculated using the Eq. ^ are depicted 
by dashed lines. It is interesting, that below a certain 
critical value k'^'' < 0.64 (corresponding to L = 2 in 
Q), the primary bifurcation is stationary indeed and the 
values of thresholds are correctly described by Eq. ([6]). 
In that case, the very first tongue crosses the p axis (in 

contrast to the situation for fco = 1.6 > fcg depicted in 
Fig. [2] where the first tongue is under the abscissa) and at 
the threshold the largest growth rate corresponds to the 
purely real branch. When approaching fcg ' from below, 

~ (c) 

this first tongue goes down and, finally, at k^ = k^ the 
whole tongue lies under the p axis. Note, that transition 
to another branch with Im\a\ ^ is accompanied by a 
stepwise change of the threshold intensity. Such a be- 
havior is related to the inclined structure of the tongues 
(see Fig. [S]). It should be stressed that for fco > 
the true values of the thresholds differs from the ones 
predicted by Eq. ^ in one order, and with increase of 
fco quite rapidly, in two orders of magnitude. The true 
threshold of the primary instability, decreases with an 
increase in fco from / = 3.3 Ic at fco = fco'^'' to 0.5 Ic at 
fco = 8 and then practically does not change at further 
increase in the fco, see Fig. |31 To obtain a stability dia- 
gram in the dimensional (L, /) plane one should rescale 
fco and p-axes of Fig. [3] multiplying them by \/{2 5n) 
and Ic, respectively. In experiment one should take ne- 
matics with rather small values of 5n to avoid very large 
values for thresholds which are of the order of Ic- For 
instance, for 5n = h ■ 10~^, K2 = 2.5 • 10~'^ dyn and 
A = 532 nm the value of Ic is 52 kW/crn^ and the critical 
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FIG. 3: Solid lines: stability diagram of the planar state on 
the (fcoiP) plane. Dashed lines on (a): thresholds calculated 
using the Eq. Solid line with fco < fco"^' (fco > ^0^^) corre- 
sponds to stationary (Hopf) bifurcation. P: region of stability 

of the plfinnr statp 




FIG. 4: Bifurcation diagrams for (a) fco = 0.5 and (b) 
fco ~ 1.6: parameter of reorientation ^' vs. p. Solid (dashed) 
curves correspond to stable (unstable) stationary solutions, 
pi: primary instability of the planar state, (a) pi — 4.76 
(stationary bifurcation); (b) pi = 0.95 (Hopf bifurcation). 
Gray region (confined by two lines): oscillatory states of the 
director. p2 = 3.65: secondary instability (heteroclinic bifur- 
cation) . 

length L'^"'^ = li^^\/{2 5n) is 34 ^m. 

It is worth noting, that the linear integro-differential 
equation ^ appears in the framework of linear analy- 
sis of another problem, namely the one with a linearly 
polarized ordinary light wave incident at a small oblique 
angle on a thin layer of homeotropically oriented nematic 
0, The only difference is that there one deals with 
the normalized incidence angle k instead of fco used here. 
It is well known that stability diagram in the (p, k) plane 
might be calculated with good accuracy in the framework 
of a two modes approximation (for k < 1) and consists of 
a line of stationary instability which joins the line of Hopf 
instability in the Takens-Bogdanov point [8] . The line of 
stationary instability exists only below a certain critical 
value of K*^"^^ and as we checked out is indeed described 
by the formula ^ after rewriting it in corresponding 
quantities. However, above k^'^^ the homeotropic state 
looses stability always via a Hopf bifurcation and the 
formula ^ is not applicable anymore. Moreover, the 
homeotropic state will never be alternatively stable and 
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FIG. 5: Dynamics near the secondary threshold p2 ~ 3.6530: 
director trajectories ai p — 3.7967 > p2 in the plane of first 
two modes ($i, $2). Starting fi-om different initial condition 
the system settles to one of the two stable nodes (filled circles). 
Red curve: separatrix. Empty circles: unstable fixed points 
(saddles). Origin: unstable spatially uniform planar state 
(focus). Dashed green line: limit cycle a.t p — 3.638 < p2- 

unstable for any k and at higher intensities a very com- 
plex dynamics occurs jH • Here we have a similar picture 
but with respect to the critical fcg . 

We choose a sum of squares of all modes ^I' = 
as a theoretical measure of reorientation inside nematic. 
It is an appropriate quantity because ^' = for the pla- 
nar state, it increases with increasing of the reorientation, 
saturates to a certain value since the amplitudes damp 
after certain n and is independent on material parame- 
ters. In Fig. m typical bifurcation diagrams are shown 
for ItQ < klf ^ and ko > A:q'^\ The planar state remains 
stable when p < pi. At p = pi we deal with a continu- 
ous transition via either a stationary for feg < A:q or a 
Hopf bifurcation in opposite case. Above the threshold 

~ (c) 

and for kg < kg the director settles to a stationary dis- 

~ (c) 

torted state, whereas for kg > k{^ to an oscillatory one. 
In latter regime, the lower and the upper lines depicted 
in Fig. Sfb) bound the region in gray and correspond to 
the minimum and maximum values taken by ^E" during 
its oscillation. The director motions develops along the 
limit cycle in the space of $„ (which is infinite). Impor- 
tantly, the planar state is always unstable for p > pi in 
both cases. In some narrow region around p2 = 3.65 the 
period of oscillations increases progressively with increas- 
ing light intensity, and oscillations become substantially 
anharmonic. The period diverges at p = p2 which cor- 



respond to a secondary bifurcation transforming the sys- 
tem into a new stationary distorted state. The dynamics 
near the threshold is summarized in the projection of the 
true phase trajectory on the plane {^1,^2), see Fig. [S] 
The bifurcation at p = p2 belongs to rather a rare type 
and corresponds to the following. Above p2 a system has 
two pairs of stationary nontrivial solutions which are mu- 
tual images under the symmetry transformation S. The 
first pair is represented by a stable node and its image 
and the second one by a saddle and its image. There 
is also a trivial solution which is represented by an un- 
stable focus corresponding to the spatially uniform pla- 
nar state. Starting from different initial conditions the 
system eventually settles to one of the stable nodes. A 
separatrix, which partitions phase plane onto the basins 
of attraction for the two stable nodes, goes through the 
two saddles and the unstable focus. Approaching p2 from 
above, the saddle and node go closer and closer to each 
other and finally merge at p = p2 . At this point the limit 
cycle appears which exist for pi < p < p2. It is worth 
noting that the symmetry S is spontaneously broken at 
the secondary bifurcation. We checked that the bifurca- 
tion scenario displayed in Fig. 3] exists over wide region 
of ko- However, the oscillatory regime seems to be more 
complex for ko > 3 than that described above. 

In conclusion, we have studied theoretically the transi- 
tions induced by linearly polarized light incident perpen- 
dicularly to a layer of nematic that has initially spatially 
uniform planar alignment. We have found the primary 
threshold as a function of (normalized) thickness of the 
layer by performing a linear analysis of the basic state. 
Our numerical analysis of the problem shows that, with 
increasing light intensity, the planar state becomes un- 
stable in a favor of either a stationary distorted state 
when the thickness is below a certain critical value, or to 
an oscillatory state if the thickness is above the critical 
value. As the intensity increases further, the oscillatory 
state disappears via a secondary bifurcation. At the bi- 
furcation point two identical saddle-node fixed points are 
born on opposite sites of the limit cycle. It results in di- 
vergence of the period of oscillations and partition of the 
limit cycle into two heteroclinic orbits. Further increase 
in the intensity brings about splitting each of the fixed 
point into a stable node and a saddle. In typical cases 
the calculated thresholds have values substantially lower 
than those believed before. 
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Miroshnichenko, E. Brasselet and A. Krekhov for helpful 
discussions and comments. 
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